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Abstract. The possibility that star-forming galaxies may leak ionizing photons is at the heart of many present-day studies 
that investigate the reionization of the Universe. We test this hypothesis on local blue compact dwarf galaxies of very high 
excitation. We assembled a sample of such galaxies by examining the spectra from Data Releases 7 and 10 of the Sloan 
Digital Sky Survey. We argue that reliable conclusions cannot be based on strong lines alone, and adopt a strategy that includes 
important weak lines such as [O i] and the high-excitation Hen and [Ariv] lines. Our analysis is based on purely observational 
diagrams and on a comparison of photoionization models with well-chosen emission-line ratio diagrams. We show that spectral 
energy distributions from current stellar population synthesis models cannot account for all the observational constraints, which 
led us to mimick several scenarios that could explain the data. These include the additional presence of hard X-rays or of shocks. 
We find that only ionization-bounded models (or models with an escape fraction of ionizing photons lower than 10%) are able 
to simultaneously explain all the observational constraints. 
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1. Introduction 

By analyzing the entire Sloan Digital Sky Survey (SDSS-III) 
DR10 spectroscopic data base (Ahn et al. 2014), we found an 
appreciable number of compact star-forming galaxies with ex¬ 
tremely high [Oiii]/[Oii] ratio^j exceeding values of 10 and 
reaching up to 50. Hereafter these galaxies are called ‘extreme 
BCDs’, standing for ‘extreme blue compact dwarf galaxies’. 

Extreme BCDs are local low-mass galaxies that are under¬ 
going intense bursts of star formation. They represent a sub¬ 
set of Green Peas (Cardamone et al. 2009) when they are in 
the redshift range 0.1 - 0.3, or of luminous compact galax¬ 
ies (Izotov et al. 2011) in a wider redshift range 0.0 - 0.6. 
Therefore, similar to Green Peas and luminous compact galax¬ 
ies, extreme BCDs share some of the properties of Lya emit¬ 
ting galaxies and Lyman break galaxies at high--;: they are com¬ 
pact, have low-mass (10 8 -10 9 M©), low-metallicity (0.1-0.2 so¬ 
lar), and high specific star-formation rates (up to 10 7 yr 1 j 
(Cardamone et al. 2009, Izotov et al. 2011). Furthermore, un- 
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1 Throughout the paper, the notations Hen, [Nii], [Oi], [On], 

[Om], [Ne iii], [Arm], and [Ariv] will stand for HenT4686, 

[N ii] T6584, [O i] T6300, [Oii]T3727, [Oiii]T5007, [Ne in] T3869, 

[Arm] T7135, and [Ariv] T4740, respectively. 


like the typical local galaxies and like Lya emitting galaxies 
and Lyman break galaxies (Nakajima et al. 2013, Shirazi et al. 
2014, Shapley et al. 2014), extreme BCDs are characterized by 
high excitation, as measured by the [Om]/[On] ratio. Thanks 
to these properties, extreme BCDs should be excellent local 
proxies for high-’ Lya emitting galaxies. 

The high [O iii]/[ 0 n] ratios can be attributed to several 
possible factors (Overzier et al. 2009, Kewley et al. 2013, 
Nakajima & Ouchi 2013, Jaskot & Oey 2013): 1) a low metal- 
licity, 2) a high ionization parameter, 3) a hard ionizing radia¬ 
tion field, 4) the presence of density-bounded H ii regions. 

This last hypothesis, if proven true, may have important 
cosmological implications in studies of the secondary ioniza¬ 
tion of the Universe, because Lyman continuum (LyC) emis¬ 
sion freely escapes the density-bounded H n regions. 

In this paper, we analyze the emission-line properties of 
local extreme BCDs with the aim of finding out whether these 
high ratios are due to significant leakage of ionizing photons. 

To construct our extreme BCD sample, we have thus 
searched the SDSS-III DR10 spectroscopic data base for star¬ 
forming emission-line objects with the highest excitation, that 
is, those with the highest [O m]/[0 n] ratios that show no sign 
of an active galactic nucleus (AGN). We have merged this sam- 
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Table 1 . Atomic data used in our abundance determinations 


Ion 

Collision strengths 

Radiative transition probabilities 

On 

Kisielius et al. (2009) 

Zeippen (1982) 

Oiii 

Aggarwal & Keenan (1999) 

Galavis et al. (1997), Storey & Zeippen (2000) 

Nii 

Tayal (2011) 

Galavis et al. (1997) 

Nem 

McLaughlin & Bell 2000 

Galavis et al. (1997) 

S ii 

Tayal & Zatsarinny (2010) 

Mendoza & Zeippen (1982) 

S hi 

Tayal & Gupta (1999) 

Froese Fischer et al. (2006) 

Cl hi 

Mendoza (1983) 

Ramsbottom et al. (2001) 

Arm 

Munoz Burgos et al. (2009) 

Munoz Burgos et al. (2009) 

Ariv 

Ramsbottom et al. (1997) 

Mendoza & Zeippen (1982) 


pie with a sample of high-excitation galaxies from SDSS DR7. 
For these extreme BCDs, in contrast to high-redshift objects, 
SDSS spectra provide a lot of information in addition to the 
intensities of the strongest lines, some of it being crucial for a 
proper interpretation of the excitation conditions. 

This paper is organized as follows: In Sect. [2] we present 
the observational database, in Sect. [3] we discuss some purely 
observational trends that lead to a zero-order interpretation, in 
Sect. [4] we present the policy we followed to interpret the ob¬ 
servational data with the help of photoionization models, and 
in Sect. [5] we compare various possible scenarios with the ob¬ 
servations. The conclusions are summarized in the last section. 


2. The observational database 

We (Y. I.) searched the SDSS/BOSS DR10 spectroscopic 
database by eye to find star-forming galaxies with high 
[Oiii]/[Oii] ratios 5). Since we were interested in star¬ 
forming galaxies, we removed objects with AGN-like spectra 
(i.e., altogether strong [O n] and [O in], broad lines, strong [O i], 
[Nil], [Sn], strong Hen, in some cases, red continuum). We 
ended up with a sample of 149 objects. All of them except one 
have [Oiii] 44363 detected. In the following this sample is re¬ 
ferred to as the DR 10 sample. 

To this sample we added the sample of star-forming galax¬ 
ies selected from SDSS DR7 (Abazajian et al. 2009) by two of 
us (Y. I. and N. G.) and described in Vale Asari et al. (2014). 
This sample is composed of all the star-forming galaxies in the 
SDSS DR7 in which the [O in] 44363 line is seen and measured 
with an accuracy of at least 25 %, and for which visual inspec¬ 
tion of the spectra and images showed that they correspond to 
genuine star-forming galaxies (i.e., excluding giant H n regions 
in spiral galaxies and active galactic nuclei). Furthermore, we 
considered only galaxies where [On]43727 is observed. With 
this criterion, all galaxies from the DR7 with z ~ 0.02 are ex¬ 
cluded. In the following this sample is referred to as the DR7 
sample. 

By combining our DR10 and DR7 samples, we obtained 
a sample of 778 objects, of which 268 have [Om]/[On] ra¬ 
tios higher than 5, and 72 have [Om]/[On] ratios higher than 
10. The remaining objects are useful for the discussion because 
they allow us to consider a wider excitation sequence and pon¬ 
der the parameters leading to the highest observed values of 
[Oiii]/[Oii]. 


The spectra from DR 10 were analyzed with the IRAF soft¬ 
ware package]^ in exactly the same way as described in Vale 
Asari et al. (2014) for the DR7 spectra. In brief, the emission¬ 
line fluxes were measured ‘manually’ using the IRAF SPLOT 
routine. Correction of hydrogen line fluxes for underlying ab¬ 
sorption and reddening correction of all lines were performed 
as described in that paper. 

In both the DR7 and DR 10 samples, the abundances of N, 
O, Ne, S, Cl, Ar, and Fe were derived from the intensities of the 
lines [On] 43727, [Oiii] 45007, [Nn] 46584, [Ne hi] 43869, 
[S ii] 46716,46731, [Sin] 49069, [Cl hi] 45518+5538, 
[Arm] 47135, [Anv] 44740, and [Fem] 44988 with re¬ 
spect to H/j. We used the temperature from [O hi] 44363/5007 
for doubly and triply ionized species and a temperature derived 
from Eq. 14 of Izotov et al. (2006) for low-ionization species. 
The ionization correction factors used to obtain elemental 
abundances from ionic abundances were those described 
in Izotov et al. (2006). The atomic data involved in the 
computations are listed in Table |T] 

The resulting abundance patterns for the DR10 and DR7 
samples are presented in Appendix [A] All the objects from 
our sample have low to moderate metallicities: 12 + log O/H 
in the range 7.5 - 8.5. This is a consequence of the fact that 
we imposed that the [O hi] 44363 line be present, implying that 
the electron temperature be sufficiently high, and thus that the 
metallicity be sufficiently low. 

3. Clues from observational trends 

In Fig. [T] we represent the position of the DR7 and DR 10 sam¬ 
ples in the classical emission-line ratio diagram [O m]/H/3 vs 
[N ii |/Ho- (the famous BPT diagram from Baldwin et al. 1991). 
Throughout the paper, DR7 objects will be represented by blue 
symbols (diamonds in this figure) and DR10 objects by red 
symbols (circles in this figure). The background gray-scale plot 
is a 2D-histogram of all the galaxies in the 7th Data Release of 
the SDSS that have a signal-to-noise ratio larger than 3 for all 
the relevant lines. The gray scale at the right of the plot indi¬ 
cates the logarithm of the number of points in each pixel. We 
see that the number of points per gray pixel thus ranges from 
1 to over 1000. The blue curve represents the empirical line 

2 IRAF is distributed by the National Optical Ashonomy 
Observatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. BPT diagram for our DR10 objects (red circles) and 
those of DR7 (blue diamonds). The background gray scale map 
represents all the SDSS DR7 galaxies with a signal-to-noise 
ratio of at least 3 in the four diagnostic lines [Oiii], H/3, [N n], 
and Her. The gray scale is in logarithm of the number of objects 
per pixel, as indicated in the bar on the right. The blue curve 
represents the Kauffmann et al. (2003) line generally used to 
distinguish star-forming galaxies from AGN hosts. 


drawn by Kauffmann et al. (2003) to distinguish star-forming 
galaxies from galaxies hosting an AGN. We see that our se¬ 
lection criteria for the DR7 and DR10 indeed samples selected 
star-forming galaxies, removing AGN hosts. We also see that 
they represent an extreme population of star-forming galaxies, 
with high values of [O m |/H/j, and [N n]/Ho' ratios that reach 
extremely low values. 

In the following, we present purely observational diagrams 
to investigate the relations between [O m ]/[0 n] and character¬ 
istic parameters of the emitting nebulae. These are the metal- 
licity, expressed in units of 12 + log O/H and the H/j equivalent 
width, EWYH/j), which measures the amount of ionizing pho¬ 
tons produced by massive stars and absorbed by the gas with 
respect to the continuum at H (5, which is mainly produced by 
the population of low-mass stars in the galaxy (which, in turn, 
is a proxy for the galaxy stellar mass). 

Other parameters worthwhile considering are L(H/3), the to¬ 
tal H/f luminosity corrected for extinction, and the density n e , 
which can be derived from the [S n] 46731/6717 line ratio. 

As is well known, two important parameters in the study 
of emission-line nebulae are the mean ionization parameter U 
(defined in Sect. 4.1 1 and the mean effective temperature of the 
ionizing radiation field, T e ff. 

In a first approximation, the abundance ratio of two ad¬ 
jacent ions of the same element, for example Ar +++ /Ar ++ , is 
given by the equation 


Ar+++ 

Ar ++ 


oc U 


Q(Ai ++ ) 
Q( H°) 


a(Ar ++ ), 


( 1 ) 


where Q(Ar ++ ) i s the number of photons above the ionization 
threshold of Ar ++ , (AH 0 ) is the number of hydrogen ionizing 


photons, and a(Ar ++ ) is the Ar ++ recombination coefficient^]. 
The Ar +++ /Ar ++ ratio is thus a function of both U and of the 
mean effective temperature r e tr of the ionizing radiation, since 


of Appendix |Q. 

The ionization parameter U cannot be measured directly 
from observations, but, assuming a value for the volume-filling 
factor e, it can be obtained from the product n e L( H/j) using Eq. 
[4] below, since L(H/3) is a direct measure of the total number 
of hydrogen ionizing photons, <2(H°) (in the hypothesis of no 
leakage and minimal absorption of LyC photons by dust). 

It has been proposed by Vflchez & Pagel (1988) that 77 , de¬ 
fined as (0 + /0 ++ )/(S + /S ++ ) is a good indicator of 7j,|) . We actu¬ 
ally to use a ratio based on line intensities rather than on ionic 
fractions, since it is readily obtained from observations. We 
also prefer to use a ratio that increases with T e ff. Thus, rather 
than the rj parameter, we use the ([O m]/[Oii])/([S iii]/[S 11]) ra¬ 
tio. As shown in the bottom right panel of Fig. |BT3] in Appendix 
|b] this ratio is far from being completely independent of U, 
and is sensitive to temperature especially at low values of T e g. 
A similar ratio is ([Ar iv]/[Ar iii])/([0 m]/[0 11]) , which is much 
less sensitive to I/, but has a lower dynamic range as a function 
of 7; ff (see the bottom left panel of Fig. |B.3| >. 

Finally, an important criterion to be considered when ex¬ 
amining the question of FyC escape is the behavior of the [0 1 ] 
line. In density-bounded objects, this line is the first to be sup¬ 
pressed, or, at least, weakened. 


Q(Ar ++ )/Q(H°) is an increasing function of 7’ c fr (see Fig. 


C.l 




Fig. 2. The [O m]/[0 11 ] ratio as a function of the oxygen abun¬ 
dance 12 + log O/H (left panel) and as a function of the equiv¬ 
alent width of the H/3 line (right panel). DR10 objects are rep¬ 
resented in red and those of DR7 in blue. The total number of 
objects in each panel is indicated in the top right corner. 


Figure [2] shows the relation between [O m]/[0 11 ] and the 
oxygen abundance (left) and EWYH/j) (right) for the DR7 and 
DR10 samples. Clearly, the highest values of [Om]/[On] cor¬ 
respond to the lowest metallicities and to the highest values of 
EW(H/I). Note that in the range of metallicities covered by our 
DR7 and DR10 samples, a decrease in metallicity from 12 + 

3 Note, however, that the effect of charge transfer reactions with 
neutral hydrogen may modulate this relation, especially at low values 
of U. 
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log O/H = 8.5 to 7.4 produces an increase in [Oiii]/[Oii] by 
about 0.3 dex if keeping the ionization parameter equal (see 
Appendix [B] Fig. |B.l| >, while the observed variation is of about 
1.3 dex. Thus, there is something special about the objects with 
the highest [Om]/[On], not only their low metallicity. The ob¬ 
served increase of [Om]/[On] with EW(H/7) suggests an in¬ 
crease with T e ff or perhaps an increase with L(Hj3), the total Hj3 
luminosity (since the continuum at H/3 is mainly due to low- 
mass stellar populations and not to the ionizing stellar pop¬ 
ulations). In any case, it argues against increasing LyC pho¬ 
ton leakage, since the latter would produce a drop in EW(H/3). 
It must be noted that at the highest metallicities encountered 
in the studied galaxies, dust absorption reduces the H/3 lumi¬ 
nosities at the highest values of the ionization parameter (see 
Appendix |B| Fig. |B.2| >. But at the low metallicities where the 
highest [O m]/[0 ii] values are found, the effects of dust absorp¬ 
tion on EW(H/7) are not important. 




Fig. 3. The [O ni]/[0 ii] ratio as a function of the electron den¬ 
sity as derived from the [S n] A6731/6717 line ratio (left panel) 
and as a function of the total luminosity in the H/3 line, cor¬ 
rected for extinction. 




log ([OIII]/[OII])/([SIII]/[SII]) log ([ArlV]/[Arill])/([OIII]/[OII]) 

Fig. 4. The [O iii]/[0 n] ratio vs ([Om]/[0 n])/([S m]/[S ii]) (left 
panel) and vs ([Ariv]/[Arm])/([0 m]/[On]) (right panel). 


trend both with increasing 7 di and increasing U in our sample 
of galaxies. The increase in 7 d r might actually be a mere 
metallicity effect, since ionizing spectra of more metal-poor 
stellar populations are harder (see Fig.|C.2|in Appendix [C|. 



Figure [3] shows the relation between [O m]/[0 ii] and two 
other fundamental properties of the emitting regions: the elec¬ 
tron density, n e , as derived from the [S n] 46731/6717 line ratio 
(left), and L(Hf3 ) (right). There is a slight tendency for high val¬ 
ues of [O m]/[0 ii] to occur at higher values of n e (the Pearson 
correlation coefficient P is 0.33). Note that the effect of colli- 
sional deexcitation at higher densities is only minor in this plot, 
since between log n e = 2 and 3.3, the emissivity ratio of [Oiii] 
and [O ii] increases by only 0.1 dex. There is no significant cor¬ 
relation between [Om]/[On] and L(H/3) {P = 0.13). 

Figure [4] shows the relation between [Om]/[On] 
and two parameters strongly related with the mean 
effective temperature of the ionizing radiation field, 
Ttff: ([Oiii]/[Oii])/([Siii]/[Sii]) (left panel) and 

([Ariv]/[Arm])/([Om]/[On]) (right panel). Both diagrams 
show a correlation, very strong for the first one (P = 0.74), 
weak but still significant for the second one CP - 0.21). 
This suggests an increase of [Om]/[On] with Z’efr- Given that 
([O m]/[0 n])/([S m]/[S ii]) is much more dependent on U than 
([Ariv]/[Ariii])/([Om]/[On]) (see Appendix[B] Fig. | B.3 [ i , Fig. 
[4] suggests that the trend observed in [O iii]/[0 n] expresses a 


Fig. 5. The [Om]/[On] ratio as a function of the product n e 
L(H/3 ) (left panel) and as a function of [O i]/[0 n] (right panel). 

Figure [5] (left) shows the relation between [Om]/[Oii] and 
the product n e L(H/3), which as shown above, can be consid¬ 
ered as a proxy for U if the filling factor e is constant. The fact 
that [O m]/[0 ii] is not seen to correlate with n e L(H/3) suggests 
that the variations in [O m]/[0 n] are due to an intrinsic increase 
of r eff rather than of U. Another option could be that e varies 
systematically from lower to higher values as [Om]/[On] in¬ 
creases. We have no way to check this directly with the data 
at hand, but this is a possibility. Figure [5] (right) shows the re¬ 
lation between [Om]/[On] and [Oi]/[Oii]. The first thing to 
note is that most objects do show [O i] emission, and that the 
[O i] emission is not limited to objects with the lowest values 
of [O m]/[0 ii] . Actually, most of the objects with the highest 
values of [Om]/[On] show [Oi]/[Oii] ratios higher than 1%. 
This suggests that the explanation for the highest values of 
[O m]/[0 ii] does not lie in LyC photon leakage. However, pho¬ 
toionization models are needed to confirm this interpretation. 
This is the subject of the next section. 
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4. Photoionization model analysis 

4.1. The photoionization model grid 

The models were constructed with the photoionization code 
Cloudy version 13.03 (Ferland et al. 2013). The abundance ra¬ 
tios with respect to oxygen are debited as follows : 

O : N : Ne : S : Cl: Ar: Fe = 1 : 0.0513 : 0.1854 : 0.0219 : 
0.00029 : 0.0048 : 0.0148, corresponding to the mean values 
of the logarithmic abundances derived for our merged DR7, 
DR10 sample (see Sect.[2]and Appendix[A]). Except for some of 
the models shown in Appendix [B] which explore the effects of 
varying O/H, we adopted for all the models a value of 12 + log 
O/H = 8, equal to the average of the measured abundances in 
our merged sample. The adopted helium abundance by number 
relative to hydrogen is 0.086, that of Mg and Si is 0.0095. 

The dust-to-gas abundance ratio is obtained by applying a 
factor 2/3 x 10 (2 21_},) to the default value in Cloudy when us¬ 
ing the ‘ism’ keyword. The expression for y is obtained from 
Table 1 (Broken law, MW case) of Remy-Ruyer et al. (2014). 
The reduction factor of 2/3 is motivated by the study of Draine 
( 2011 ). 

For most models, the ionizing radiation held is given by 
the population synthesis code PopStar (Molla et al. 2009) for 
a Chabrier (2003) stellar initial mass function and at the ap¬ 
propriate metallicity, obtained by interpolation. The reference 
case corresponds to an age of 1 Myr. Other ionizing radiation 
belds have also been considered, as described in the following 
sections. All the models assume a constant hydrogen density, 
but two geometries are considered. They are debited by two 
different values of fs = R ln /Rs. where R m is the model inner 
radius and Rs is the Stromgren radius corresponding to R m = 0. 
Models with f s = 0.03 are equivalent to a full sphere, models 
with fs = 3 correspond to a hollow spherical bubble. The mod¬ 
els are computed until the ratio of ionized hydrogen to total 
hydrogen density falls below 0.02. 

In addition to the spectral energy distribution (SED) of the 
ionizing radiation, the gas chemical composition and density, 
one needs to debne the luminosity of the ionizing source. It is 
more convenient, actually, to bx the volume-averaged ioniza¬ 
tion parameter, U estimated in the following way. In a nebula, 
the ionization parameter U(R) at a distance R from the center 
is debned as 


U(R) = 


Q( H°) 
4nR 2 nc ’ 


( 2 ) 


where (AH 0 ) is the ionizing luminosity in photons s , n is the 
hydrogen density and c the speed of light. The volume averaged 
ionization parameter, U, is debned by 


_ f U(R)dV 
U = ~r - 

fdV 


(3) 


where the integration is performed over the volume of the ion¬ 
ized nebula. The radius of the ionized nebula is obtained by 
equating Q(H°) with the total number of hydrogen recombina¬ 
tions to upper levels in the ionized volume. Assuming that the 
gas is uniformly distributed in the nebula with a constant den¬ 


sity n and a volume hlling factor e, and assuming that the case 
B recombination coefficient, o- B , is constant, we obtain 

d input = -jr ( Q(H°)ne 2 ) ((1 + /|) 1/3 - fs). (4) 

Our series of models assume a density n = 100 cm~ 3 and 
take values of log f/i nput ranging from -4 to -1 by steps of 0.5 
dex, which allows us to span the entire range of observed values 
of [Oiii]/[OiiQ As we show below, slight differences occur in 
the predicted line ratios with the same value of U i npu t for the 
two geometries considered. 

For each model, we created several density-bounded ver¬ 
sions by trimming the outer parts and leaving only the inner 
parts corresponding to / H /j = 20, 40, 60, 80, and 100% of the 
value of the H/j luminosity obtained for the radiation-bounded 
model. 

4.2. How the models are used 

In the present work we aim to determine the cause of the high¬ 
est values of [O m]/[0 n] observed and, in particular, whether 
they can be attributed to LyC photon leakage. One of the im¬ 
portant lines to consider is therefore [O i] 46300, which is pro¬ 
duced in the warm transition region between fully ionized gas 
and neutral gas. The diagnostic diagram we chose to consider 
is [Oiii]/[Oii] versus [O i]/[0 in]. In a brst approximation, it is 
independent of O/H, which allows us to compare the models 
with 12 + log O/H = 8 with our entire observational sample. 

Since we suspect that the highest [Om]/[On] might be 
linked to high values of T e g, it is of interest to use lines of ions 
with high ionization potentials as diagnostics. One such line is 
that of [Ariv], since the ionization potential of Ar ++ is 40.6 eV, 
signibcantly higher than the ionization potential of 0 + (35.1 
eV). Another option would be the [Ne m] line, since the ioniza¬ 
tion potential of Ne + is also 40.6 e V, and, in addition, there is no 
charge transfer between Ne ++ and Ne + , in contrast to what oc¬ 
curs between Ar +++ and Ar ++ . The advantage of [Ariv] is that 
it can be related to [Arm], and that the [Ariv]/[Ar hi] ratio does 
not depend on abundances. The second diagnostic diagram we 
therefore chose to consider is [O m]/[0 n] versus [Ar iv]/[Ar m]. 

If we wish to understand the behavior of emission-line 
ratios in galaxies, we cannot put aside some lines, such as 
Hen, which, although weak, provide important diagnostics be¬ 
cause they indicate a very hard ionizing radiation field (above 
54.4 eV). Their systematic presence in Hii galaxies has first 
been noted and commented onby Garnett et al. (1991) and 
Stasinska & Izotov (2003), and was more thoroughly discussed 
by Shirazi & Brinchmann (2012) and Jaskot & Oey (2013). 
All the suggested explanations (Wolf-Rayet stars, shocks, X- 
ray binaries, stellar rotation effects on the effective temperature 
of massive stars, inaccurate predictions of stellar atmosphere 
models at high energies) can be disputed. In the present pa¬ 
per, our aim is not to propose an astrophysically valid expla¬ 
nation, so we consider different scenarios to explain the ob¬ 
served He n/H/j ratios and examine their effects on the remain- 

4 In the insets of the figures we replaced log(/i nput by log U for 
simplicity. 
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ing line ratios. We therefore choose as a third diagnostic di¬ 
agram [O m]/[0 n] versus He n/H/3, and retain photoionization 
models that predict He n/H/j- 0.01, which is representative of 
the bulk of the objects in our sample that show Hen emission 
(the He n line is not detected in about one third of the objects of 
our sample, but in at least part of these cases this absence can 
be attributed to a too low signal-to-noise ratio). 


5. Comparing models with observations 

5.1. Models using classical stellar population 
synthesis results 

5.1.1. The reference model: PopStar with an age of 1 
Myr 

Figure [6] shows the distribution of the observational 
points in our three diagnostic diagrams: [O iii]/[0 n] versus 
[Ar iv]/[Ar hi] (left), [Om]/[Oii] versus [Oi]/[Oiii] (middle), 
and [Oiii]/[Oii] versus He n/H/i (right). Superimposed is our 
grid of photoionization models computed with the SED from 
PopStar models at an age of 1 Myr and with an oxygen abun¬ 
dance 12 + log O/H = 8. The dashed lines link models with 
the same values of t/i nput , while the continuous lines link mod¬ 
els with the same value of the trimming parameter / H /j. The 
upper row of panels correspond to our filled sphere models 
(Js = 0.03), while the lower row corresponds to our bubble 
models (f$ = 3). It can be seen that these models fall very 
short of reproducing the observed He n/H/3 ratios, even if con¬ 
sidering density-bounded models. The SED from the PopStar 
models at an age of 1 Myr simply does not provide enough 
He + ionizing photons. The observational points fall well on 
the ionization-bounded models line in the [O m]/[0 n] versus 
[O i]/[0 hi] diagranj^] both for the filled and hollow sphere 
models. However, this cannot be considered as proving that the 
objects are all ionization-bounded, since the models fail to re¬ 
produce the He n/H/i and the [Ar iv]/[Ar hi] line ratios (espe¬ 
cially for the hollow sphere models). We note that the highest 
observed [O in]/[0 n] ratios correspond to the highest values of 
U input considered in our models. We also note that the hollow 
bubble models seem to provide a poorer fit to the observations 
in the [O m]/[0 ii] versus [Ar iv]/[Ar m] diagram than the filled 
sphere models, but this question cannot be fully addressed until 
the SED problem is solved (see below). 


5.1.2. Other stellar populations 

In Fig. C.3 of Appendix |C| we show the SEDs obtained at an 
age of 1 Myr with the same stellar initial mass function, but 
using the different assumptions for the stellar evolution mod¬ 
els proposed by starburst99 (Leitherer at al. 1999): Geneva 
tracks with standard mass loss, Geneva tracks with high mass 
loss, original Padova tracks, and Padova tracks with AGB stars. 


5 It is interesting to note that, in the [O m]/[0 ii] versus [O i]/[0 m] 
plot, the curves for g < 1 almost overlap, with only the /h^ = 1 curve 
standing out. This is because in that case the [O i] emission is domi¬ 
nated by the ionization front, while in the density-bounded cases it is 
produced by the traces of neutral hydrogen in the whole H n region. 


All the SEDs have been obtained with the stellar atmosphere 
grid of Smith et al. (2002) based on atmosphere models of 
Pauldrach et al. (2001) and Hillier & Miller (1998), similarly to 
the PopStar models. From this figure, it is clear that at an age 
of 1 Myr, all these SEDs are equivalent and will lead to almost 
identical photoionization models. 

Another question to consider is the age of the ionizing stel¬ 
lar populations. Fig. |C.4| of Appendix |C| shows the SEDs ob¬ 
tained by PopStar for several ages. We see that, for an age 
of 2 Myr, the SED is almost identical to that at 1 Myr, and 
that for 3 and 5 Myr, the SEDs are softer. This implies that 
they will lead to the same problems regarding [Ar iv]/[Ar iii] 
and He n/H ft. However, at an age of 4 Myr, which is when mas¬ 
sive Wolf-Rayet stars begin to come into play, the SED is much 
harder. 

Figure [7] shows the results of photoionization models com¬ 
puted with PopStar for an age of 4 Myr, using the same presen¬ 
tation as Fig. [6] We see that the [Om]/[On] versus [Oi]/[Oiii] 
diagram is again compatible with ionization-bounded models, 
but that the SED at this age is even too hard to account for the 
bulk of the observed He n/H/j ratios. One then might think that 
the relevant ages for reproducing the He n/H/1 ratios is larger 
than 3 Myr and somewhere in the vicinity of 4 Myr. 

However, this cannot be the explanation for all the objects 
that show He n emission. The He u line has been observed in 
488 objects out of the 778 objects in our sample. In a first ap¬ 
proximation, in absence of leakage or dust absorption of ion¬ 
izing photons, and assuming that all the photons with energies 
higher than 54.4 eV are absorbed by He + (which is roughly 
correct except at low values of the ionization parameter, see 
Stasiriska & Tylenda 1986), we have 

HeII/H/3 = AQ(He + )/Q(H°), (5) 

where A is the ratio of emissivities of the He n and H/( lines 
divided by the ratio of the case B recombination coefficients 
o'P,(He )/o'B(H°). The coefficient A depends only weakly on 
the electron temperature and is equal to 1.74 at 20,000 K. The 
left panel of Fig. [8] shows the values of g(He + )/Q(H°) pre¬ 
dicted by PopStar models of metallicities Z = 0.004,0.008 
and 0.02 (which correspond to 12 + log O/H =7.7, 8 and 8.4) 
as a function of age. We see that Q(He + )/Q(H°) leads to a de¬ 
tectable Hen only during an age range of 3.9—4.25, 3.4-5.4, 
and 4.0-5.25 Myr, respectively, for the three metallicities con¬ 
sidered, if we assume detectability at He ii/H/j > 10 5 . The 
right panel of Fig. [8] shows the values of He ii/H/j as a func¬ 
tion of EW(HyS) for the same PopStar models, taking Eq. [5] 
to estimate He n/H/3. This figure can be directly compared to 
the observations, shown in the left panel of Fig. [9] having in 
mind the observed distribution of metallicities as a function of 
EW(H/j), shown in the right panel of Fig. [9] We see that the 
models do not predict any Hen emission for EW(H/j) above 
200 A, while many of our objects (and especially some of those 
with the highest values of [O hi]/[0 n] (see Fig. [2]) show higher 
values of EW(H/j). We must note, in addition, that the com¬ 
puted values of EMTH/i) correspond to an instantaneous star- 
burst, while in real objects, part of the continuum comes from 
old stellar populations that do not contribute to the ionization, 
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log [ArlV]/[ArlII] 







Fig. 6. The [O m]/[0 n] versus [Ar iv]/[Ar hi] (left), [O in]/[0 n] versus [O i]/[0 in] (middle) and [O m]/[0 u] versus He n/H/i (right) 
diagnostic diagrams. The red and blue points represent the DR10 and DR7 objects respectively. The curves represent our grid of 
models computed with the SED from PopStar models at an age of 1 Myr and with a metallicity 12 + log O/H = 8. The dashed 
lines link models with the same values of t/i npu t, the continuous lines link models with the same value of the trimming parameter 
fu/ 3 , the ionization-bounded models corresponding to / H /; = 1. The key for the colors is given in the bottom right panel. The 
upper row of panels corresponds to our filled sphere models (fs = 0.03), while the lower row corresponds to our bubble models 
(fs = 3). 






Fig. 7. Same as Fig. [6] but for models at an age of 4 Myr. 


so that the observed values of EWiH/j) are lower than the the- even at values of EW(H/7) lower than 200 A, it is likely that in 
oretical values for the most recent starburst. This means that, many of the objects the observed He n line is not due to ioniza- 
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Age (Myr) EW(Hfl) 


Fig.8 . Left: (7(He f )/Q(H !l ) as a function of age for PopStar 
SEDs with metallicities Z =0.004, 0.008, and 0.02 as indicated 
in the inset. Right: log He n/H/f versus EWtH/f), where He ii/H/j 
is taken equal to 1.74g(He + /Q(H°), for the same PopStar 
models. 


tion by a simple stellar population. Considering the possibility 
of extended starbursts rather than instantaneous ones does not 
help since it dilutes the He ii emission simultaneously with the 
H/i emission. A similar discussion can be found in Shirazi & 
Brinchman (2012) who, in addition, found that 40 per cent of 
the galaxies in their sample do not have Wolf-Rayet features in 
their spectra despite showing strong nebular He n emission. 

The above discussion of the production of nebular Hen 
emission by stellar populations of about 4 Myr relies on 
PopStar models. Fig. |C. 6 | of Appendix |C| shows the SEDs ob¬ 
tained at 4 Myr for the same stellar synthesis models as Fig. 


C.3 We see that, at this age, the SEDs differ significantly owing 


to the different stellar evolutionary tracks involved. However, it 
can be intuited that the above discussion remains qualitatively 
valid and that the He n emission observed in the objects of our 
sample cannot be in majority caused by Wolf-Rayet stars. 

In the next subsections, we explore various scenarios based 
on modified ionizing SEDs that might better explain the ob¬ 
served trends in our diagnostic diagrams. 


5.2. Models using blackbodies as proxies for a harder 
stellar radiation field 

Here we consider simple blackbodies as proxies for the SEDs 
of the ionizing stellar populations. We find that a black- 
body with r e ff=60,000 K leads to a value of ~ 0.01 for 
He ii/H/3 in ionization-bounded models. Figure [10] shows the 
results of models built with this ionizing radiation field for 
our three diagnostic diagrams (i.e., it is the analog of Fig. 
[ 6 ] for this SED). We see that the values of [Oi]/[Oiii] fall 
on the ionization-bounded line, but that, again, the predicted 
[Ar iv]/[Ar hi] line ratios are too low. Considering slightly 
cooler blackbodies and density-bounded models, while accept¬ 
able for the [O m]/[0 n] vs He n/H/i diagram would be incom¬ 
patible with the [O iii]/[0 ii] vs [Ar iv]/[Ar hi] diagram, since for 
a given value of [O hi]/[0 ii] the models lead to lower values of 
[Ar iv]/[Ar hi] as T e 3 - decreases, as shown in Fig. El 

We remark that the [Om]/[On] vs [Ar iv]/[Ar 111 ] diagram 
demands a harder SED than that of a blackbody at 60,000 
K, especially at the highest values of [O iii]/[0 11] . From Fig. 




Fig. 9. Left: observed values of log Hen/H/3 versus EW(H/3) 
(for objects with no He 11 detection log He 11 /H/? has been set to 
-4); Right: observed values of 12 + log O/H versus EW(H/i). 


11 one can estimate that the SED above 40.8 eV should pro¬ 


vide roughly as many photons as a blackbody at 90,000 K 
or more. Current synthetic stellar population models do not 
consider stellar atmospheres with such a hard radiation held. 
However, it is to be noted that some of the models for very 
massive stars at low metallicities computed with an appropri¬ 
ate treatment of radiation-driven winds (Kudritzki 2002) do 
provide the amount of photons needed to reproduce both the 
[Oiii]/[Oii] vs [Ariv]/[Arm] and the [Om]/[On] vs He ii/H/l 
diagrams. Perhaps using more realistic model atmospheres to 
construct SEDs from stellar population models could help. 




Fig. 11. The [Om]/[On] vs [Ariv]/[Arni] diagram for 
ionization-bounded photoionization models with a blackbody 
ionization source (12 + log O/H = 8 ). Top panel: models 
for filled spheres; bottom panel: models for spherical bubbles. 
Continuous lines correspond to different values of the effective 
temperature, while dashed lines correspond to different values 
of the ionization parameter. The color coding is indicated in the 
plot. The positions of the observational points for our samples 
(red for DR 10, blue for DR7) are shown for comparison. This 
diagram shows that for a given value of [O iii]/[0 11 ] the value 
of [Ar iv]/[Ar hi] increases with T e ff. 
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Fig. 10. Same as Fig. [6] but with models computed with a blackbody ionizing source at a temperature of 60,000 K. 






Fig. 12. Same as Fig. [fi] but with models in which a bremsstrahlung at 10 8 K is added to the stellar ionizing radiation from 
the PopStar models corresponding to an age of 1 Myr, in such a proportion that 1% of the ionizing photons come from the 
bremsstrahlung. 


5.3. Models using classical stellar populations and an 
additional bremsstrahlung 

Another possibility to explain the observed values of He u/H/3 
might be the presence of massive X-ray binaries (Van Bever 
& Vanbeveren 2007) or of shocks in the atmospheres of stars 
leading to production of X-rays (Pauldrach et al. 2001), or else 


X-ray emission from protostars (Bonito et al. 2010). We mim¬ 
icked all these possible scenarios by adding a bremsstrahlung 
component with a temperature of 10 8 K to the PopStar mod¬ 
els of Sect. 5.1.1 This temperature is typical of the values in¬ 
ferred from a detailed analysis of A 0535+26, one of the best- 
studied massive X-ray binary (Rothschild et al 2013). We ad- 
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justed the contribution from the bremsstrahlung so as to ob¬ 
tain He n/H/j- 0.01 in the resulting photoionization models. 
This implies that 1% of the ionizing photons come from this 
bremsstrahlung component of the spectrum. The resulting SED 
is shown in Fig. C.5 of Appendix [C] The results of the pho¬ 
toionization models in our diagnostic diagrams are shown in 
Fig. 12 (to be compared to Fig. [6). We see that, at the high¬ 
est values of U the ionization models slightly deviate from 
the observations in the [O m]/[0 n] vs [O i]/[0 hi] plot, but that 
density-bounded models reproducing the observations have a 
LyC photon escape fraction of 10% or lower. We note, how¬ 
ever, that in the [Oiu]/[On] vs [Ariv]/[Arm] plot the obser¬ 
vational points tend to lie below the models, meaning that the 
ionizing source does not provide enough photons at energies 
slightly above 40.6 eV. 


5.4. Models using classical stellar populations and 
adding the effects of shocks 


Another option to consider is the effects of shocks either pro¬ 
duced by supemovae from previous generations of stars or from 
random motions of interstellar clouds. Such shocks could eas¬ 
ily produce He n lines and were even advocated to explain the 
[Ne v] /13426 emission observed in some star-forming galaxies 
as predicted by photoionization models (Izotov et al. 2012). We 
adopted the strategy of Izotov et al. (2012) and Jaskot & Oey 
(2013) to add the contribution of shocks to the emission-line 
properties of star-forming galaxies. We used the grid of radia¬ 
tive shock models from Allen et al. (2008) and extracted from 
it the models constructed with the SMC abundance set because 
this chemical composition is the most appropriate for our case 
among the available compositions. 

Figure [13] shows the variations of [Ar ivj/H/j, [O rJ/H/k and 
He n/H/j as a function of the shock velocity for three different 
values of the magnetic field. The photoionization models and 
shock models are combined in the following way: 


ic/m = 


(Iph/Hfiph + I s h/H(3 sh x C) 
(1+C) 


( 6 ) 


where I/H[3 stands for the intensity of a line relative to H/T the 
subscripts ‘c’, ‘ph\ and ‘sh’ stand for ‘combined’, ‘photoion¬ 
ization’, and ‘shock’, respectively, and C is the adopted ratio of 
L{H/3) sh /L{H/3) ph . 

In Fig. [14] we show the combination of the models using 
PopStar with a shock model with a velocity of 300 km s 1 and 
a magnetic field strength of 0.5 /jG. This value of the mag¬ 
netic field minimizes the contribution to He n/H/3 and maxi¬ 
mizes the contribution of [O r|/H/f The value of C to obtain 
He n/H/i= 0.01 is 0.09. This figure shows that, at the highest 
values of [O hi]/[0 ii] the combined models produce an [O i \/H/3 
ratio significantly higher than observed, regardless of the value 
of the trimming parameter fiyj. The reason is that the [O i] emis¬ 
sion in a combined model that explains the He n/H/3 is com¬ 
pletely dominated by the shocks. At low values of [O iii]/[0 n], 
on the other hand, the effect of the shocks on [O m]/[0 n] is 
imperceptible. 

In Fig. [15] we show the combination of the photoioniza¬ 
tion models using PopStar with a shock model with a veloc¬ 


ity of 300 km s _I and a magnetic field strength of 10 pG, this 
time maximizing the contribution to He n/H/3 and minimizing 
the contribution to [O i]/H J3. The value of C to obtain He n/H/3 
= 0.01 is 0.055. In this case, the grid of models follows the 
observational points in the [O m]/[0 n] versus [O i]/[0 hi] plot 
rather well, with however a slight deviation from the observed 
trend at the highest values of [O m]/[0 n]. 

We thus find that for some combinations, a fraction of the 
observations are in fact compatible with density-bounded mod¬ 
els, implying a leakage of ionizing photons in these cases (but 
the [Ar iv]/[Ar hi] ratio remains unsatisfactory). Note that at the 
highest values of [Om]/[On] the [Oi]/[Oiii] ratio is entirely 
determined by the shock and that the curves corresponding to 
different values of /h/j merge. Estimating the amount of FyC 
photon leakage would at least require a simultaneous fitting of 
all the relevant line ratios for each object and an independent 
estimate of the ionization parameter. This procedure would in¬ 
volve many unconstrained parameters, such as the shock veloc¬ 
ities or the strength of the magnetic field. 

If, however, unlike what is assumed in our combined mod¬ 
els, the shocks occur in the ionized region, the [O i] emission 
would be different. From basic knowledge of the physics of 
ionized nebulae, one can intuit that the extreme UV and X-ray 
photons produced by the shock would travel until the outskirts 
of the H ii region and produce an [O i] emission similar to that 
of our PopStar + bremsstrahlung models. 

From these considerations, we conclude that the pres¬ 
ence of shocks can enhance the He n/H/3 ratio and cause it to 
agree with the observations and, at the same time, enhance 
the [O i]/[0 m] ratio, especially for the objects of highest ex¬ 
citation. In these circumstances, an analysis based only on 
emission-line ratios cannot determine the existence or impor¬ 
tance of FyC photon leakage. 

IWe note, however, that in the case of the objects with high¬ 
est excitation the [O i]/[0 hi] ratio is essentially determined by 
the shock. If the ionized nebulae in these objects were density- 
bounded, the tight correlation observed between [Om]/[On] 
and [O i]/[0 hi] would imply a very strong dependence of the 
shock properties on the characteristics of the photoionized neb¬ 
ula, in particular its ionization parameter. This is not impossi¬ 
ble but is an additional requirement for the hypothesis of FyC 
photons leakage to be the correct interpretation. 

5 . 5 . Models including the effect of an active nucleus 

The final scenario we consider is that of an AGN embedded 
in a star-forming region. Again, the SED of the stellar popula¬ 
tion is provided by the PopStar models at an age of 1 Myr. For 
the AGN we took the standard AGN SED provided by Cloudy. 
This AGN spectrum was combined with the PopStar spectrum 
so as to obtain He n/H/3 - 0.01. The contribution of the AGN to 
the total number of ionizing photons must be quite small oth¬ 
erwise the He n/H/3 ratio is boosted toward values much higher 
than observed (this contribution is about 3% to obtain He n/H/3 
= 0.01). The combined SED is shown in Fig. |C.5| of Appendix 
C] The results of the photoionization models are shown in Fig. 
16] Here, the models clearly indicate that the observed nebulae 
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Fig. 13. The line ratios He ii/H/i (left panel), [O i l/H/i (middle panel), and [Ar iv|/H/i (right panel) as a function of shock velocity 
for the shock models from Allen et al. (2008) (see text). The different curves correspond to different values of the magnetic field, 
with the color code indicated in the left panel. 






log [Oll/lOIII] 


log Hell/H/3 



log [OI]/[OIII] 



Fig. 14. Same as Fig. [h] but with combined photoionization and shock models. The photoionization models use the energy 
distribution of PopStar at an age of 1 Myr. The shock models are extracted from the shock-model library of Allen et al. (2008) 
for a magnetic field of 0.5 /iG and have a velocity of 300 km s 1 . They are added in such a proportion to obtain He wjWji- 0.01 
for the combined model (see text for more details). 


would be ionization-bounded. As in Figs. [6] and [12] the models presence of photons at high energies (40.8 eV and 54.4 eV, re¬ 
give [Ar iv]/[Ar hi] ratios somewhat lower than observed. spectively) that might blur the information from the [Oi] lines. 


6. Summary and discussion 

We have tested the hypothesis of LyC photon leakage from lo¬ 
cal extreme BCDs by assembling a sample of high-excitation 
galaxies from the SDSS DR 10 and DR7 and by examining 
their spectroscopic properties in the light of photoionization 
models. We have argued that a proper diagnostic cannot rely on 
strong lines alone: some weak lines are of prime importance to 
provide an unambiguous answer. Such are the [Oi] line, emit¬ 
ted in the warm transition region between ionized and neutral 
gas, and the [Ar iv] and He n lines, which are indicative of the 


After showing that photoionization models using SEDs 
from currently available stellar population synthesis models - 
even at ages when Wolf-Rayet stars are present - cannot ac¬ 
count for all the observational constraints for the entire sam¬ 
ple, we considered several scenarios that could explain the data. 
These include the additional presence of a hard X-ray spectrum 
that might be produced by massive X-ray binaries or shocks in 
the atmospheres of stars or protostars or by an active galac¬ 
tic nucleus. Another option would be shocks that might be 
produced by supernovae or cloud-cloud collisions. It turns out 
that only ionization-bounded models (or models with an escape 
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Fig. 15. Same as Fig. Il4l but for shock models with a magnetic field of 10 /./G. 






Fig. 16. Same as Fig. [fi] but with models where a AGN-type ionizing spectrum is added to the stellar ionizing radiation from 
the PopStar models corresponding to an age of 1 Myr. The contribution of the AGN radiation is chosen to produce a resulting 
He n/FT/3 of =* 0.01. 


fraction of Lyc photons lower than 10%) are able to explain all 
the observational constraints at the same time. 

We also found that the current stellar population synthe¬ 
sis models provide too soft ionizing radiation fields to explain 
the observed [Ar iv]/[Ar hi] ratios. This is probably because the 
stellar atmosphere models incorporated in the codes underpre¬ 


dict the number of photons emitted between 40.8 and 54.4 eV, 
and perhaps even at higher energies. 

Our study also showed that the extreme BCDs have higher 
ionization parameters (reaching log U- -1.5 or even -1 de¬ 
pending on the nebular geometry) than galaxies with lower val¬ 
ues of [O iii]/[ 0 n] ratios. Such high values of the ionization pa- 
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rameters are probably produced by a more compact distribution 
of the gas around the ionizing stars. 

The question of the possible leakage of LyC photons in 
high-redshift galaxies should be examined with methods simi¬ 
lar to tthose presented in this paper, and not relying on strong 
line ratios alone. This will be possible with future deep sur¬ 
veys such as the Keck Baryonic Structure Survey (Steidel et al. 
2014). 

We must acknowledge, however, that our study is based on 
idealized photoionization models. Perhaps some leakage might 
occur because of a nebular covering factor smaller than one and 
not because the nebulae are density-bounded? Perhaps the [O i] 
emission actually comes from optically thick clumps embed¬ 
ded in a density-bounded medium? Such hypotheses need to 
be tested by taking into account additional considerations (for 
example, imaging studies of nearby galaxies). In the meantime, 
our Occam’s razor approach of photoionization modeling, the 
fact that [Oi]/[Oiii] correlates so strongly with [Om]/[On], 
and the fact that the highest [Om]/[On] ratios correlate with 
the highest H/i equivalent widths argue against strong LyC 
photon leakage from the galaxies with the highest observed 
[Oiii]/[Oii] ratios. 
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Appendix A: Chemical abundances 

Figure [AT] shows the abundance ratios N/O, Ne/O, S/O, Ar/O, 
Cl/O, Fe/O as a function of O/H for both our DR 10 and DR7 
samples (red and blue points, respectively). Note that we mea¬ 
sured the S abundance only for objects from our DR10 sam¬ 
ple since in the DR7 sample, the [S hi] 49069 line for all the 
objects with measured [On] 43727 lies beyond the spectro¬ 
scopic range. The slight observed tendency between S/O and 
O/H might be due to inaccurate ionization correction factors. 
An improvement of the ionization correction factors presented 
in Izotov et al. (2006) would require a dedicated study that is 
beyond the scope of the present paper. The large dispersion in 
the computed Cl/O abundance ratios is essentially due to the 
extreme weakness of the lines and the large associated error 
bars. The dispersion in the N/O ratios is real and discussed in 
detail in Vale Asari et al. (2014). The trend of Fe/O decreasing 
with increasing O/H, already noted by Izotov et al. (2006) us¬ 
ing galaxies from SDSS DR3, is an indication that more Fe is 
depleted into dust grains as metallicity increases. 

For the photoionization modeling described in Section[4]we 
took a chemical composition defined by the mean abundances 
of our DR7-DR10 sample. 

Appendix B: Didactic diagrams 
Appendix C: Spectral energy distributions 
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Fig.B.l. The effect of metallicity on the [O iii]/[0 ii] ratio ob¬ 
tained from photoionization models. The models are com¬ 
puted for a PopStar ionizing radiation field (age 1 Myr) and 
are ionization-bounded. Continuous lines join models of same 
metallicity, while dashed lines join models with same ioniza¬ 
tion parameter. The color code is given in the plot. The posi¬ 
tions of the observational points for our samples (red for DR 10, 
blue for DR7) are shown for comparison. 





12 + log 0/H 



Fig.A.l. Abundance ratios in our galaxy sample. DR10 objects 
are plotted in red, DR7 objects are in blue. The total number of 
objects included in each panel is indicated in the upper right 



logU 


Fig. B.2. The values of L(H/5)/(9(H 0 ) as a function of log C/ input 
as computed in ionization-bounded photoionization models 
with different metallicities. The ionizing radiation field is from 
PopStar at an age of 1 Myr. Continuous lines join models of 
same metallicity, while dashed lines join models with same 
ionization parameter. The color code is given in the plot. This 
diagram shows how the absorption of the ionizing radiation by 
dust grains affects the H/J luminosity depending on the dust 
abundance (which is linked to the metallicity, see Sect. |4.1| > 
and on the ionization parameter. 
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Fig. B.3. The variations of several line ratios with black- 
body temperature for different values of the ionization param¬ 
eter (ionization-bounded models with 12 + log O/H = 8). 
Continuous lines join models of same effective temperature, 
while dashed lines join models with same ionization parame¬ 
ter. The color coding is the same as in Fig. 11 




Fig. C.2. Same as Fig. |C. 1 1 for SEDs obtained by PopStar at 


an age of 1 Myr for a Chabrier (2003) stellar initial mass and 
different metallicities Z as indicated by the color code in the 
inset. 



C.l 


for several SEDs obtained with star- 


Fig. C.3. Same as Fig. 
burst99 with a Kroupa (2001) stellar initial mass function for 
a metallicity Z = 0.004 at an age of 1 Myr. The color coding 
indicates the different cases proposed by starburst99: Geneva 
tracks with standard mass loss, Geneva tracks with high mass 
loss, original Padova tracks, and Padova tracks with AGB stars. 
The PopStar model at an age of 1 Myr is shown for reference. 
We see that, at this age, all the SEDs are exactly sumperim- 
posed. 


Fig. C.l. The log of Q(E), the number of photons emitted 
above the energy threshold £ as a function of E, for different 
spectral energy distributions (SEDs). The ionization thresholds 
of important ions are indicated by the vertical dashed lines. The 
different curves correspond to blackbodies at different temper¬ 
atures, as indicated in the inset. The PopStar model at an age 
of 1 Myr of Fig. C.4 is shown for comparison. 
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Fig. C.4. S ame as Fig. |C.l| but with SEDs corresponding to Fig. C.6. Same as Fig. |C.3| but for an age of 4 Myr. Here, all 


PopStar models at a metallicity Z = 0.004 for different ages 
(1,2, 3, 4, and 5 Myr). 


the SEDs are harder than in the 1 Myr case, because of Wolf- 
Rayet stars. We also note a substantial difference between the 
various SEDs. 



E [eV] 


Fig. C.5. S ame as Fig. |C.4| for SEDs corresponding to AGN 


spectral energy distribution, combined PopStar at an age of 1 
Myr and AGN (Sect. [53), and combined PopStar at an age of 
1 Myr and bremsstrahlung (Sect. 5.3 i. The pure PopStar model 
at an age of 1 Myr of Fig. C.4 is shown for comparison. 
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